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I.   INTRODUCTION  AND  BACKGROUND 

Early  Developments 

Commencing  three  years  ago,  a  feasibility  study  was  undertaken  to  determine 
if  an  effective  aircraft  fatigue  monitoring  instrument  could  be  conceived  on 
the  basis  of  microprocessors.   The  goal  was  to  design  a  small,  lightweight 
package,  which  could  record  maximum  and  minimum  loads  sequentially  in  a  form 
conducive  to  automated  data  processing.   The  first  unit  was  designed  and 
fabricated  by  LT  David  Vidrine  (1)  and  LT  James  Sturgis  (2).   This  first  pro- 
totype was  physically  larger  than  would  have  resulted  for  an  optimum  design 
solely  for  fatigue  monitoring  in  order  to  give  it  testing,  check-out,  and 
development  capability.   It  was  built  around  the  Intel  8008  CPU  chip,  and  had 
16  channel  capacity  via  a  multiplexer.   It  digitized  the  incoming  load  infor- 
mation from  strain  gages,  accelerometers,  or  other  transducers  and  then 
searched  these  signals  for  maxima  and  mimima.   Upon  finding  them,  the  digital 
value  of  each  extremum  point  was  stored  in  memory  in  chronological  order  until 
the  memory  was  filled,  at  which  point,  the  microprocessor  caused  the  tape 
recorder  to  be  activated,  and  the  memory  to  be  transferred  to  the  tape  cassette. 

This  unit  was  laboratory  tested  with  sinusoidal  signals  from  signal 
generators  with  results  as  shown  in  Fig.  1.   After  success  was  obtained  with 
known  signals,  the  unit  was  flight  tested  in  an  NPS  leased  Cessna  310.   Typi- 
cal results  are  shown  in  Figure  2. 
Second  Prototype 

In  those  budding  days  of  microprocessors,  the  state  of  the  art  was 
advancing  rapidly  as  it  is  today,  and  by  the  time  the  first  prototype  was 
built,  new  developments  had  been  perfected.   LT  Wesley  Stanfield  (3)  began 
work  on  a  second  prototype  that  was  directed  toward  flight  test  in  a  high 


performance  aircraft.   It  was  based  on  the  8008  CPU,  but  incorporated  ad- 
vancements in  programming  capability  and  reduction  in  the  physical  size  of  the 
components.   It  was  designed  with  8  channel  capacity  and  had  the  physical 
dimensions  shown  in  Figure  3.   The  final  version  of  the  software  was  generated 
by  LT  Charles  Butler  (4) ,  who  also  generated  the  programs  for  data  reduction 
from  the  instrument  through  a  Hewlitt  Packard  9830,  desk-top  computer. 

This  second  prototye  was  successfully  laboratory  tested  and  demonstrated 
for  Dr.  David  Weiss  of  the  Naval  Air  Development  Center  (NADC) .   The  instrument 
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was  able  to  take  data  from  signal  generators  up  to  5H   ,  pick  out  maxima  and 
minima  and  record  them  on  tape.   Using  the  HP  9830,  the  voltage  of  the  maxima 
and  minima  are  listed  in  tabular  form  in  chronological  sequence. 
Final  Results 

This  concluded  the  feasibility  study  at  NPS,  and  the  Naval  Air  Systems 
Command,  through  Bob  Weinberger,  adopted  the  concept  of  an  in-flight  micro- 
processor instrument  to  monitor  fatigue  as  well  as  other  flight  and  engine 
parameters.   The  F-18  was  designated  as  the  first  Navy  aircraft  to  incorporate 
the  concept,  and  all  future  aircraft  will  use  state-of-the-art  versions  of  this 
fatigue  monitoring  concept. 

The  second  prototype  instrument  was  then  to  be  transitioned  to  Dr.  Weiss 
at  NADC  for  further  experimentation  while  new  prototypes  were  designed  and 
packaged  to  meet  military  specifications.   A  malfunction  in  the  analog  to 
digital  conversion  unit  was  experienced  at  this  time.   Steps  taken  to  correct 
it  finally  required  a  new  unit  from  an  alternate  manufacturer,  which  in  turn 
precipitated  changes  in  the  wiring  and  associated  changes  in  the  software.   In 
the  meantime,  further  developments  in  microprocessor  technology  have  made 
significant  investments  of  time  in  the  second  prototype  unit  unwarranted. 


Future  Possibilities 

In  the  original  prototyping,  a  tape  recorder  was  used  as  the  final 
reservoir  of  the  data.   This  was  a  temporary  arrangement  because  of  the  limited 
conditions  under  which  successful  tape  recorder  operations  can  occur.   Large 
temperature  and  humidity  variations  can  cause  the  capstans  to  freeze,  and 
recording  in  a  high  g  environment  may  not  be  reliable  at  all;  however,  there 
is  an  alternative  emerging. 

Two  specific  possibilities  have  recently  come  forth:  one  is  a  CMOS  storage 
cell,  which  requires  so  little  power  that  small  batteries  can  sustain  them 
into  the  time  ranges  of  weeks  that  would  be  needed  for  actual  operations . 
The  second  is  Texas  Instrument's  bubble  memory,  which  at  this  time  is  not 
completely  reliable  but  which  has  promise  to  fill  the  needs  of  the  system  in 
the  not  too  distant  future. by  providing  large  quantities  of  nonvolatile  memory 
capacity. 


II.   DAMAGE  CALCULATION 

In  light  of  the  Navy  commitment  to  employ  microprocessor  units  as  fatigue 
monitors,  it  now  becomes  appropriate  to  look  forward  to  the  type  of  data  that 
will  be  available  from  them  and  how  it  can  be  used  to  calculate  fatigue 
damage. 
Two  General  Approaches 

The  Navy  for  many  years  has  used  a  cumulative  damage  approach  to  fatigue, 
based  upon  Miner's  Law  or  some  variation  of  it,  which  has  now  been  developed 
into  a  rather  complex  calculation  that  considers  various  methods  of  counting 
the  cycles,  plastic  flow  around  the  point  of  stress  concentration,  residual 
stresses  and  their  relaxation,  etc.   On  the  other  hand  the  Air  Force  is  pur- 
suing fatigue  calculations  based  upon  crack  propagation.   An  initial  flaw  of 
.050  in.  is  as summed  to  exist  at  each  critical  stress  riser  location,  and  the 
number  of  cycles  required  to  grow  it  to  a  critical  size  constitutes  the 
fatigue  life. 

Two  large  computer  programs  have  been  imported  and  made  operational  on 
the  NPS  IBM  360.   Potter's  (5)  cumulative  damage  program  uses  Miner's  Law  and 
is  presently  operational.   Engels  (6)  program,  "Cracks  II,"  employ  the  crack 
propagation  method  and  is  now  being  checked  out  with  example  problems  provided 
by  the  author.   Check  runs  are  planned  in  the  future  to  give  comparisons  of 
fatigue  lives  for  the  same  conditions  and  to  investigate  the  affect  of  cycle 
minimums  in  the  load  spectrum.   The  calculations  will  be  for  an  actual  air- 
plane fatigue  critical  point  on  the  A-7E. 
Potter  Program 

LT  Scott  Atkinson  (7)  has  investigated  the  influence  of  counting  method 
and  type  of  block  loading  on  damage.  The  load  spectrum  used  in  the  analysis 
was   the  spectrum  which  is  reproduced  in  Table  I,  from  MIL  SPEC  A8866,  and  the 


TABLE  I 

FREQUENCY  OF  MANEUVER  LOADS 

NUMBER  OF  TIMES  PER  THOUSAND  HOURS  THAT 
LOAD  FACTOR  IS  EXPERIENCED 


PERCENT  OF  MAXIMUM  (POSITIVE)                FLIGHT  MANEUVER  LOAD 

SYMMETRICAL  LIMIT  LOAD  FACTOR  SPECTRUM  A 

35  17000 

45  9500 

55  6500 

65  4500 

75  2500 

85  1360 

95  440 

105  150 

115  40 

125  16 

TOTAL  42006 


PERCENT  OF  MAXIMUM  (NEGATIVE)                FLIGHT  MANEUVER  LOAD 

SYMMETRICAL  LIMIT  LOAD  FACTOR  SPECTRUM  A 

0  500 

10  200 

20  100 

30  60 

40  35 

50  30 

60  25 

70  20 

80  15 

90  10 

100  5 

110  3 

TOTAL   1003 


loads  are  assumed  to  be  applied  at  wing  station  32  on  an  A-7E  wing  (8).   Load 
data  were  input  as  block  loads  and  also  as  single  randomized  cycles.   The 
Miner's  Law  summation  of  fatigue  damage  for  this  loading  spectrum  is  .1300 
per  1000  flight  hours.   Results  obtained  using  this  computer  program  are 
given  in  Tables  II,  and  III.   It  is  important  to  notice  that  when  fatigue 
damage  is  calculated  without  using  the  range-pair  counting  method,  it  makes 
only  a  minor  difference  on  the  total  fatigue  damage  when  the  sequence  of 
loading  is  changed  or  when  negative  loads  are  included. 

Block  loads  were  arranged  in  L0-HI,  HI-LO,  and  HI-LO-HI  sequences.   The 
results  indicate  the  L0-HI  sequence  to  be  more  damaging  than  the  HI-LO 
sequence,  with  damage  due  to  the  HI-LO-HI  sequence  falling  between  the  other 
two  block  load  sequences.   In  the  HI-LO  sequence  the  HI  loads  leave  a  local 
residual  compressive  stress  in  the  material,  and  the  L0  loads  do  less  damage. 

In  the  case  of  randomly  generated  load  sequences,  identical  starter 
integers  were  used  in  order  to  evaluate  the  effect  on  the  damage  of  inclusion 
of  MIL  spectrum  A  negative  loads  and  ground  cycles.   The  more  negative  loads 
that  are  included,  the  worse  is  the  fatigue  damage.   The  increase  in  damage 
with  more  negative  loads  is  caused  by  the  negative  loads  decreasing  the  local 
residual  compressive  stress  caused  by  positive  loads  into  the  plastic  zone. 
This  permits  the  local  positive  elastic  stresses  to  be  more  damaging. 
Minimum  Load  Levels 

Including  the  negative  values  of  stress  resulting  from  landing  loads, 
which  through  the  range-pair  counting  method  is  tantamount  to  appropriately 
accounting  for  the  ground-air-ground  cycle,  makes  considerable  difference  in 
damage  accumulated  even  for  this  high  performance  airplane.   This  is,  of  course 
contrary  to  many  statements  in  the  literature  that  GAG  cycles  are  only  of 
significance  in  transport  aircraft,  and  this  conclusion  is  proper  when 


TABLE  II 


SEQUENCE  OF 
BLOCK  LOADS 

LO-HI 

HI-LO 

HI-LO-HI 


FATIGUE  DAMAGE  DUE  TO  POSITIVE  MILSPEC  A 
BLOCK  LOADING,  NORMALIZED  TO  1000  HOURS 


11%  LIMIT  LOAD  MIN. 
USING  RPCM* 

.16780102 

.069243014 

.10864216 

*   RANGE-PAIR  COUNTING  METHOD 


11%  LIMIT  LOAD  MIN 
NOT  USING  RPCM* 

.16780108 

.069243014 

.10839599 


TABLE  III 

FATIGUE  DAMAGE  DUE  TO  MILSPEC  A 
RANDOMIZED  LOAD  INPUTS,  USING  RPCM,  NORMALIZED  TO  1000  HOURS 


INTEGERS  USED 
TO  START  RANDU 


83745,54711,54487 

13547,66549,7 

9,583,4777 

48621,3,491 

73,559,1001 

357,833,1 


MINIMUM 

LOAD  VALUES 

11  PERCENT 

NEGATIVE 

GROUND 

LIMIT  LOAD 

MILSPEC  A 

CYCLES 

.039226338 

.046267733 

.060848035 

.095124505 

.085063167 

.10429338 

.043321513 

.043640956 

.036781987 

.018672124 

.034077277 

.035881512 

.031819174 

.042006299 

.037130679 

.037665060 

.049336702 

.050973855 

residual  stresses  are  not  being  accounted  for.   However  the  negative  loads 
cause  a  reduction  in  the  protective  residual  stress  at  the  stress  concentration 
site  and  an  associated  increase  in  damage  as  shown  by  the  third  column  of 
Table  III. 

The  entire  question  of  minimum  loads  is  raised  at  this  point.   Since 
minimum  loads  are  not  recorded  at  each  cycle  on  the  accelerometer  counters, 
it  has  always  been  assumed  that  the  minimum  of  each  cycle  was  lg. ,  which  is 
intuitively  thought  to  be  unrealistic.   However  from  the  data  in  Table  III, 
the  minimums  are  important  and  that  information  will  be  recorded  and  available 
for  use  from  the  new  microprocessor  fatigue  monitoring  instruments.   In  the 
meantime,  flight  testing  on  current  aircraft  with  the  same  instruments  can 
be  used  to  generate  a  spectrum  of  minimum  load  levels  for  each  cycle.   This 
spectrum  can  then  be  randomized  and  used  in  the  program  in  the  same  way  that 
the  11%  loads  and  the  negative  MILSPEC  A  loads  were  used  in  this  study. 
Random  Load  Histories 

Figure  4  illustrates  the  pattern  of  the  individual  load  randomization 
process  and  portrays  the  first  two  flight  hours  of  the  aircraft.   The  range  of 
values  in  the  columns  of  Table  III  are  caused  by  such  variations  in  the  load 
randomization  pattern. 

The  fatigue  damage  resulting  from  randomization  of  loads  on  a  cycle  basis, 
as  shown  in  Table  III,  is  significantly  different  from  the  damage  due  to  the 
HI-LO-HI  block  loads  illustrated  in  Table  II,  which  has  been  accepted  by  many 
as  being  representative  of  random  loading.   For  instance,  Vought  Aeronautics 
Division  recently  conducted  a  fatigue  study  using  randomized  flights  where 
each  flight  was  based  on  a  mission  profile  specified  by  the  Air  Force.   They 
found  fatigue  damage  resulting  from  the  randomized  flights  was  in  close  agree- 
ment with  fatigue  damage  resulting  from  HI-LO-Hi  block  loading. 


Table  IV  has  been  compiled  to  show  the  influence  of  not  using  the  range- 
pair  counting  method,   van  Dijk's  (9)  arguments  are  so  compelling  concerning 
the  advisability  of  using  the  range-pair  method  and  these  results  show  that  it 
should  be  used  although  it  is  considerably  more  expensive  to  run.   The  range- 
pair  counting  method  relates  the  load  sequence  to  the  stable,  cyclic,  stress- 
strain  behavior  of  a  material  by  coupling  up  strain  ranges  so  that  each 
combination  of  half  cycles  will  produce  a  complete,  closed  stress-strain  hyste- 
resis loop  in  order  to  be  counted  as  a  cycle. 

The  range-pair  counting  method  is  illustrated  in  Figure  5.   It  counts  a 
strain  range  as  a  cycle  if  it  can  be  paired  with  a  subsequent  straining  of 
equal  magnitude  in  the  opposite  direction.   For  a  complicated  load  history, 
some  of  the  ranges  counted  as  cycles  will  be  simple  ranges,  such  as  2-3, 
during  which  the  strain  does  not  change  direction,  but  others,  such  as  1-8, 
will  be  interrupted  by  smaller  ranges,  which  will  also  be  counted  as  cycles. 
In  Figure  5,  ranges  are  marked  with  solid  lines  and  the  paired  ranges  with 
dashed  lines. 

Each  peak  is  taken  in  order  as  the  initial  peak  of  a  range,  except  that 
a  peak  is  skipped  if  the  part  of  the  history  immediately  following  it  has 
already  been  paired  with  a  previously  counted  range.   If  the  initial  peak  of 
a  range  is  a  minimum,  a  cycle  is  counted  between  this  minimum  and  the  most 
positive  maximum  which  occurs  before  the  strain  becomes  more  negative  than 
the  initial  peak  of  the  range.   For  example,  in  Figure  5  a  cycle  is  counted 
between  peak  1  and  peak  8,  peak  8  being  the  most  positive  maximum  that  occurs 
before  the  strain  becomes  more  negative  than  peak  1.   If  the  initial  peak  of 
a  range  is  a  maximum,  a  cycle  is  counted  between  this  maximum  and  the  most 
negative  minimum,  which  occurs  before  the  strain  becomes  more  positive  than 
the  initial  peak  of  the  range.   For  example,  in  Figure  5  a  cycle  is  counted 


between  peak  2  and  peak  3,  peak  3  being  the  most  negative  minimum  before  the 
strain  becomes  more  positive  than  peak  2.   Each  range  that  is  counted  is 
paired  with  the  next  straining  of  equal  magnitude  in  the  opposite  direction, 
explaining  why  complete  cycle  rather  than  half  cycle  counts  are  made.   For 
example,  in  Figure  5  part  of  the  range  between  peaks  8  and  9  is  paired  with 
the  range  counted  between  peaks  1  and  8. 
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INTEGERS  USED 
TO  START  RANDU 


TABLE  IV 

FATIGUE  DAMAGE  DUE  TO  MILSPEC  A  RANDOMIZED  LOAD 
INPUTS,  NOT  USING  RPCM,  NORMALIZED  TO  1000  HOURS 

MINIMUM  LOAD  VALUES 


11  PERCENT 

NEGATIVE 

GROUND 

LIMIT  LOAD 

MILSPEC  A 

CYCLES 

83745,54711, 

54487 

.040602274 

.047398917 

.049124956 

13547,66549, 

7 

.040137805 

.047727041 

.049612150 

9,583,4777 

.041204430 

.048346408 

.049540661 

48621,3,491 

.041160211 

.049230121 

.050957575 

73,559,1001 

.040792227 

.047043338 

.049311966 

357,833,1 

.041056350 

.049220584 

.050481893 
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III.   LOCAL  STRESS  CALCULATION 

There  are  two  aspects  of  this  program  that  require  further  investigation, 
and  the  first  of  these,  local  stress  calculation  will  be  discussed  in  this 
section.   By  local  stress  is  meant  the  stress  at  the  site  of  the  stress 
concentration.   The  transfer  function  between  the  local  stress  and  the  far-fielc 
strain  of  a  skin  in  uniaxial  tension  must  be  a  function  of  the  geometry  of  the 
notch  or  hole,  the  local  strain,  and  the  type  of  material. 


a  =  f  (E,p,  a,  e,  s,  e)  (1) 


where 


E  =  Young's  Modulus 

p  =  Geometric  parameter  -  radius  of  curvature 
of  the  notch  or  hole. 

a  =  Geometric  parameter  -  size  of  the  notch  or 
hole 

e  =  Local  strain 

s  =  Far-field  stress 

e  =  Far-field  strain 

To  nondimensionalize  this  set  of  variables,  stress  concentration  factors 
are  used.   For  instance  any  notch  can  be  fit  with  an  ellipse  as  proposed  by 
Inglis  (11) ,  and  the  stress  concentration  factor  for  an  ellipse  in  the  elastic 
range  is 


Kt  «  7  -  [1  +  2  /  -]  (2) 


Defining  dimensionless  ratios  for  local  stress  and  strain  in  the  plastic 
range  in  the  natural  way, 

K   5  -    k  =   ~  (3) 

as     e   e  w/ 
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the  transfer  function  can  be  simplified  to 

K0  -  /  (K£,  KT)  (4) 

Neuber  (12)  proposed  a  function  to  describe  this  relationship  in  the  form 
of  a  geometric  mean  of  the  stress  and  strain  concentration  factors. 

K  K  =  K,,,2  (5) 

a  e    T 

This  was  evaluated  on  7075-T651  experimentally  by  Garske  (13)  by  using  large 
notched  plates  so  that  strain  gages  could  be  placed  at  the  stress  concentration 
for  reliable  local  strain  readings.   Four  plates  were  studied  with  K  =  3.92  , 
3.33  ,  2.51  and  1.63  .   Neuber' s  relationship  was  found  to  be  in  error  in 
predicting  local  stresses  and  strains  by  as  much  as  24%  in  the  higher  plastic 
range  in  all  four  sheets. 

From  the  appearance  of  the  data  in  graphical  plots  and  guided  by  Neuber 's 
original  proposal,  a  slightly  generalized  form  of  the  power  law  was  proposed. 

K0X  K/  =  A  KT2  (6) 

To  properly  maintain  the  elastic  limit  condition  where  K  =K  ■»  K  ,  A  ■  1 
and  z  =  x  +  y  . 

K/  K/  -  **  (7) 

Without  loss  of  generality,   y  can  be  set  to  one  for  later  ease  in  calculation, 
Rearranging, 

fejfe)'1 
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A  curve  fit  of  the  stress  and  strain  concentration  in  the  plastic  range  was 

made  of  the  data  for  plate  #3  for  which  Km  =  2.51.   For  a  correlation 

v  T 

coefficient  of   .90,    the  result  was 

/K   X1'244   /K    \ 

~     -    -966  (9) 


*T/ 


T> 


The  fact  that  this  came  very  close  to  the  proper  elastic  limit  constant  of 
unity  was  very  encouraging  that  the  general  form  is  correct.   Whether  there  is 
something  fundamental  in  the  form  of  the  equation  and  1.24  should  be  5/4  is  not 
known. 

It  has  the  general  appearance  of  an  energy  density,  but  after  several 
attempts  to  derive  it  from  fundamental  principles  failed,  it  was  taken  as  an 
empirical  law  as  it  is.   However  to  keep  the  appropriate  and  consistent  elastic 
form,  the  equation  should  be 

/K  X1*244  /K  \ 

te)    te)  ■ 1 

In  further  calculations  this  equation  (10)  will  be  used. 

To  carry  out  the  calculations,  the  stress-strain  relationship  was  formulate 

from  experimental  daca  on  7075-T651  in  the  classical  Ramburg-Osgood  form  by 

a  power  law  curve  fit  of  the  plastic  stress-strain  data  points. 

f     1  22.678 
+  99§15  (11) 


Ll0.22xl06J 

This  expression  fit  the  nonlinear  portion  of  the  uniaxial  data  with  a  correla- 
tion coefficient  of  .98. 

The  final  expression  of  local  stress  as  a  function  of  far-field  strain  can 
now  be  formulated  from  the  rearranged  geometric  equation  (10) . 


14 


f) 1-244  (t)  -  s2-2" 


Since     s  =   Ee   , 


(12) 


(I) 


1.244 


(i^e) 


2.244 


(13) 


Inserting  the  stress-strain  expression  for  e  , 


,2.244 


ylO.  22x10 


a. 27x10 


n       \23.922         ... 
°—)  "  (KTe)2'244  -  0 


(14) 


Solutions  of  this  equation  for  stress  over  a  range  of  values  of  far-field 
strains  gave  maximum  variations  of  4%  from  the  test  data  for  the  three 
different  K  ' s   tested. 

With  the  local  stress  determined  accurately,  the  residual  stress  can  be 
determined  accurately  by  assuming  elastic  unloading.   For  instance,  if  the 
fatigue  monitor  returned  with  a  sequence  of  recorded  strains  as  shown  in 
Table  1,  the  corresponding  local  stresses  could  be  calculated  by  the  techniques 
that  have  been  developed  here. 


TABLE  I 

KT 

=  1.63 

Cycle 

Far-field 
strain 

Far-field 
stress  (psi) 

Local  Stress 
(psi) 

1st  Peak 

.005308 

54,248 

77,163 

1st  Valley 

.000000 

0 

-  11,260 

2nd  Peak 

.001455 

14,870 

12,960 

2nd  Valley 

.000000 

0 

-  11,260 
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The  first  peak  local  stress  comes  from  equation  (14)  using  the  given  far-field 
strain  and  K  .   The  far-field  stress  producing  it  is  s  =  E  (.005308)  =  54,2 41 
The  first  valley  stress  is  calculated  by  removing  the  far-field  stress  of 
54,248  psi,  which  corresponds  to  locally  removing  Aa  =  K  S  =  (1.63) (54,248)  = 
88,424  psi.   Subtracting  88,424  psi  from  the  peak  stress  of  77,163  psi  gives 
the  residual  stress  of  -  11,260  psi  after  the  first  cycle.   The  second  cycle 
strain  of  .001455  produces  a  stress  increase  from  equation  (14)  of  24,226  psi. 
However  since  there  was  a  residual  of  -  11,260  psi,  the  second  cycle  peak 
stress  will  be  12,966  psi.   Upon  unloading,  since  there  has  been  no  plastic 
yielding,  the  structure  returns  to  the  same  residual  stress  state  as  before. 
The  first  cycle  yielding  has  been  of  considerable  benefit  in  protecting  the 
structure  during  the  second  cycle  with  a  compressive  residual  stress,  which 
reduces  the  fatigue  damage  accumulated.   This  can  be  repeated  throughout  a 
randomized  or  an  actually  recorded  load  history  where  the  range-pair  counting 
method  has  been  incorporated. 
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IV.  RESIDUAL  STRESS  RELAXATION 

The  second  of  the  two  aspects  of  the  program  that  requires  investigation  is 
the  relaxation  of  the  residual  stress  created  by  loads  into  the  plastic.  In  the 
true  sense  of  the  word,  stress  relaxation  is  stress  decay  under  constant  strain, 
but  for  the  present  the  term  stress  relaxation  will  be  used  more  loosely  to  mean 
stress  reduction  with  time  and  will  be  descriptive  of  residual  stress  behavior 
with  time. 

When  structures  containing  stress  concentrations  are  loaded  sufficiently  to 
produce  local  plastic  straining  in  the  region  of  high  stress,  residual  compres- 
sive stresses  remain  in  that  region  upon  unloading  of  the  structure,  as  has  been 
described.   The  residual  stress  and  strain  are  highly  localized  disturbances  in 
the  elastic  continuum,  and  since  they  are  disturbances,  they  are  potentially 
unstable.   Given  sufficient  conditions,  they  will  tend  to  be  reduced  and  will 
relax  toward  lower  values  during  succeeding  cycles  until  the  next  overload  occurs. 
Potter  (10)  has  proposed  that  the  residual  stress  be  decomposed  into  an  equilib- 
rium portion  and  a  transient  portion  and  suggested  an  exponential  to  model  the 
transient  behavior,  thinking  of  it  as  being  similar  to  a  transient  in  a  critically 
damped  dynamic  system.   The  overload  triggers  the  impulse  and  controls  its  height, 
and  following  the  impulse,  there  is  an  asymptotic  return  to  an  equilibrium  state. 
Potter's  model  for  residual  stress,   aD  ,  could  be  expressed  as  follows: 

aR  =  aREQ  +  aRTRANS  (15) 

In  terms  of  the  proposed  exponential, 

bN 
aR  "  °REQ  +  Ae  (16) 

N  =  Number  of  cycles 

where  it  is  assumed  that  the  transient  residual  stress  decreases  with  cycles 
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A  and  b  are  evaluated  from  the  initial  conditions  and  the  end  conditions.   If 


a     =  o         at  N  =  0,   then 


A  ■  °R0  "  °REQ  (17) 


If  N  =  N    is  the  point  when  the  transient  residual  stress  is  90%  relaxed;  i.e. 
EQ 


aR  =  aR  "  -9  (aR0  "  °REQ> 


b  =  %^i  (18) 

^EQ 


this  gives  for  a  , 
R 


,  .  .  .In  0.1  N 

°R  =  aREQ  +  (°R0  "  °REQ)e(       ^  >  (19) 


or 


aR  =  aRE0  +  (0RO  "  aREQ)e  "  ^  V  (20) 

Experimental  tests  on  uniaxial  specimens  were  run  by  Bentley  (14)  to  eval- 
uate the  undetermined  parameters  in  this  model  and  study  this  general  behavior. 
The  first  area  of  interest  was  in  the  mechanism  of  the  substructural  recovery,  b? 
relaxation  of  the  lattice  structure,  that  accounts  for  the  stress  decay  toward 
a  stabilized  stress  level.   Cyclic  relaxation  tests  on  7075-T7651  into  the  plast: 
zone  between  fixed  strain  levels  revealed  no  detectable  stress  relaxation.   Con- 
sistently every  relaxation  test  showed  the  reloading  stress— strain  curve  rejoin* 
at  the  point  of  departure  from  the  original  stress— strain  curve  for  the  cyclic 
unloading,  which  is  evidence  of  zero  relaxation.   However,  another  recovery 
mechanism,  anelastic  behavior,  was  measured,  and  it  was  concluded  that  it  was 
the  primary  element  producing  residual  stress  decay. 

Anelastic  behavior  is  illustrated  in  Figure  6,  where  it  is  seen  that  all  of 
the  plastic  strain  present  at  unloading  is  not  permanent  set,  but  some  of  it  is 
recovered.   After  the  load  was  removed,  the   £    continued  to  decrease  exponen- 


tially  with  time  for  periods  of  the  order  of  one  hour.   Using  the  same  form  of 
the  model  as  we  did  for  stress, 

E  =  £    • ,  +  e„    •  (21) 

equxl.    transient 

and 

e^aT,a  =  ^n  ~   e  )exp(-2.303  t/t  )  (22) 

trans     p    eq  eq 

At  a  point  of  stress  concentration  in  the  structure,  the  material  under 
highest  stress  may  be  thought  of  as  a  miniature  uniaxial  specimen  as  shown  in 
Figure  7.   After  plastic  straining  at  this  location,  the  specimen  is  longer 
than  it  was  originally,  and  the  surrounding  elastic  material  is  trying  to  com- 
press it  back  to  its  original  length.   This  puts  the  element  in  compression  and 
is  where  the  residual  stress  arises  from.   However,  as  the  plastic  strain 
decreases  through  anelastic  recovery  the  mismatch  is  gradually  reduced  and  not 
as  great  as  at  first.   The  surrounding  material,  consequently,  does  not  exert 
as  great  a  compressive  stress  as  it  did  originally.   The  transient  portion  is 
given  by 

*!,«.     =  (aDn-aD   )exp(-2.303  t/t   )  (23) 

Rtrans     RO   Req  eq 

This  process  was  found  by  Bentley  to  be  time  dependent  rather  than  cycle 
dependent.   The  anelastic  recovery  was  measured  on  one  specimen  undergoing  a 
single  load  cycle,  which  was  then  left  to  sit.   It  was  then  measured  on  a  specimen 
which  was  cycled  continuously.   The  anelastic  recovery  was  the  same  in  each  case. 
The  residual  stress  decay  is  a  function  of  time  and  not  cycles. 

The  total  expression  for  residual  stress  is 

CR  ■  °Re<,  +  (°R0  -  V)eXP('2J°3  t/teq)  (24) 

If  a  is  expressed  as  a  percentage  of  a        ,  the  anelastic  recovery  tests 

can  be  used  to  evaluate  the  unknown  parameters  in  the  model.   The  o*    is  pro- 

RO 
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duced  by  the  mismatch  in  the  strains  at  the  edge  of  the  notch 


aR0  =  k  £P 


(25) 


The  equilibrium  stress,   a    ,  is  produced  by  the  mismatch  in  strains  remaining 

eq 

after  time. 


o_   =  k  e   =  cr_n  ( — 3.) 
Req      eq    RO   e 

P 


(26) 


Using  these  in  the  expression  for  a_  yields 

K 


r  ,  e 


aR  =  aR0 


11 


+     1 


_  -Bl 


exp  -  2.303   t/t 


eq 


(27) 


Tests  were  run  at  two  different  strain  levels  in  the  plastic  zone  with  results 
as  given  in  the  following  table: 

TABLE  II 


P 

e 
eq 

e  U 
eq  p 

t 
eq 

.00540 

.005033 

.932 

53.49  min 

.005046 

.004784 

.948 

50.05  min 

A  broader  range  of  strain  levels  in  the  plastic  zone  would  need  to  be  run  to  deter 

mine   if   e  Iz        is  a  constant  and  if   t    is  a  material  constant, 
eq   p  eq 

If  this  were  the  case,  the  residual  stress  expression  would  become 


°t>   =  a™  [°-94  ~  °-06  exP  ("  0.0919t)] 
K     RO 


where  o_   is  found  as  described  in  the  last  section. 

KO 


(28) 


Certainly  these  few  tests  cannot  be  relied  upon  to  accept  this  expression 
as  factual  for  7075-T651;  however  two  important  points  have  been  forthcoming. 
First,  the  stress  decay  was  dependent  upon  time  and  not  upon  cycles,  and  secondly 
it  was  small  in  magnitude.   All  of  this  will  be  probed  more  thoroughly  and  deeply 
during  the  next  year. 
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V.   SUMMARY  AND  CONCLUSIONS 

A  fatigue  monitoring  instrument  based  upon  a  microcomputer  has  been  demon- 
strated in  laboratory  and  preliminary  flight  test  to  be  feasible  and  very  desir- 
able.  Due  to  its  adoption  for  future  aircraft  by  NAVAIR,  investigations  have  been 
begun  for  effective  use  of  the  far-field  strain  data  recorded  by  it  to  compute 
fatigue  damage. 

Three  major  areas  have  been  undertaken  in  this  study:  (1)  A  sophisticated 
computer  program  has  been  used  to  investigate  the  effects  of  the  minimum  portion 
of  each  load  cycle  on  the  accumulated  fatigue  damage.  (2)  Analytical  and  expe- 
rimental investigations  have  been  made  of  the  use  of  Neuber's  equation  in  calcu- 
lating local  stress  from  far-field  strain  with  the  result  that  a  better  equation 
was  developed.  (3)  An  experimental  investigation  of  the  relaxation  of  residual 
stresses  produced  at  sites  of  stress  concentration  was  undertaken  with  some 
success  but  with  other  work  to  be  done  before  these  results  are  conclusive. 
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Figure  1  Laboratory  Test  Results  of  the 
Fatigue  Monitor 
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Figure 


3  MIDAS  Aircraft  Fatigue  Monitor 
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Figure  5      Illustration  of  Range-Pair   Counting  Method 
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Figure  6   An  Illustration  of  Anelastic  Behavior 
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Figure  7  Uniaxial  Representation  of  Residual  Stress 
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